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Abstract 

We present a method for symbolic calculation of Feynman amplitudes for 
processes involving both massless and massive fermions. With this approach 
fermion strings in a specific amplitude can be easily evaluated and expressed 
as basic Lorentz scalars. The new approach renders the symbolic calculation of 
some complicated physical processes more feasible and easier, especially with the 
assistance of algebra manipulating code for computer. 
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1 Introduction 



In high energy colhsions, generally, multi-particle final states are produced. To calculate 
a practical physical process with several fermions in the final state, conventionally, is a 
tedious work to implement, even at the tree level. One needs to square the Feynman 
amplitude, sum over fermion polarizations, and to take trace for each possible fermion 
string loops. The non-Abelian nature of the Standard Model makes the number of 
independent Feynman diagrams grow very rapidly as the external particles increas- 
ing. Fortunately, due to the recent development in computer technology, a lot of time 
consuming and tedious Dirac algebra and symbolic simplification work can be carried 
out on a ordinary personal computer. Nevertheless, even with the up-to-date algebra 
manipulating programs for computer, many a multi-particle process tend to be still 
cumbersome to calculate. On the other hand, another drawback of conventional ampli- 
tude squaring technique is that it loses the information on spin correlations, which is 
more and more attainable in present experimental measurements. 

An alternative technique, the so-called helicity method, has been developed pretty 
applicable in practice in recent years. With this method, the independent fermion 
strings are contracted and reduced to analytic expressions of four-vector products, such 
that the matrix-element squaring becomes a trivial work. The development of helicity 
method has experienced a bit long history. Some early efforts on this can be seen in 
Refs. 12 121 E]- In recent development, a bunch of approaches for calculating the helicity 
amplitudes are proposed in Refs. [4-14], and it should be noted the references given 
here are far from complete on this subject. 

In practice, most of physical processes involve fermion production and decays, which 
give a number of independent fermion lines within a Feynman amplitude. The generic 
form of a fermion line involving a pair of fermions can be expressed as 



where (pi,Ai) and (^2)^2) stand for the momenta and polarizations of the external 
fermions; and in between of the two spinors there is a string of Dirac matrices. Since real 
physical processes are described by the matrix element square, in conventional method 
of calculating an interaction process, one squares the amplitude, which contains more 
or less (QJ-like fermion lines when there are external fermions involved, and sums over 
the fermion polarizations to get fermion loop(s). Using the nature of spinor projection 
operators, one then is confronted with couple of traces. As aforementioned, with the 
growth of the number of external fermions one need to evaluate a vast number of fermion 
loop traces. Even with the help of computer algebra manipulation, quite often one finds 
it is still a cumbersome work to perform the calculation. 

Hehcity method enables one to evaluate physical amphtude directly, resulting in 
a limited number of four-vector scalars. In literatures, there are several kinds of ap- 
proaches in evaluating the fermion strings analytically. One of them is carried out by 
choosing a specific representation for the Dirac matrices [0], and then get the explicit 
forms of all the related matrices. Another kind of popular helicity method fOJ is 
performed by introducing extra lightlike auxiliary momenta and expressing the massive 
particle spinors by the corresponding massless ones. 





In this work, we use the maximum credits of Refs. pTlj ITT] to generahze a straight- 
forward method for the evaluation of fermion strings. We beheve this approach, which 
leads to a covariant form for the helicity amplitude, is a relatively simple and fast way 
in practical computations. 



2 Formalism 



In the conventional method one computes unpolarized cross sections by squaring the 
eq. ([Q)-like amplitude while summing over fermion polarizations, and obtains 

tr[(^i ± mi) fS'i ^2 ■ ■ ■ ^„(j^2 ± ma) ^„ ■ ■ ■ ^2 ^1] , (2) 

where mi and m2 are the masses of pi and p2, and the plus and minus signs before 
them correspond to particle and anti-particle, respectively. Taking the above trace in 
principle is not a tedious work, especially with the help of computer. However, with 
the number of external fermions increasing, in conventional amplitude modulus square 
approach one encounters a prohibitive number of eq.Q-like traces to evaluate, which 
makes the physical computations very time consuming, even for the machine. 

Noticing that the fermion string (Q) can be rewritten as a trace form, like 

tr[^i ^2 • ■ ■ iSMP2, A2) ® u{Pi: Ai)] , (3) 

therefore, now, the question of how to evaluate the helicity amplitude becomes how to 
re-express u{p2,X2) ® ^{Pi,Xi) by basic Dirac matrices. In Ref.[13j, Vega and Wudka 
obtained one kind of expression for it by making use of the Bouchiat-Michel identity 
While in principle their result is applicable for practical calculation, we still feel 
it not so convenient for programming, since the main result of (9) in Ref.jTB] is not 
expressed in a convenient and covariant way for different helicities. 

The helicity method we will present is a straightforward generalization of the results 
of Refs.pni El- By choosing two auxiliary four- vectors, /cq and ki, one can fix up the 
basic spinors u{ko, A) by the equations: 

u{ko, X)u{ko, X) = ^ ^^'^^ h (4) 

and 



u{kQ, X) = X }iiu{kQ, —A) . (5) 

Here, 

ko ■ ko = 0, ki ■ ki = —1, A;o ■ /ci = . (6) 

With the above defined basic spinors, the wavefunction of a massive fermion (particle 
or anti-particle) can be expressed as 

^± m 

nv, X) = ^ ^ u{ko, -A) . (7) 

^/2p ■ ko 



It is easy to check that these spinors, and their conjugations defined by the normal way, 
satisfy Dirac equations; they are also the eigenstates of 75 fi with eigenvalues of ±1. 
Here, n is the fermion polarization vector. Using Eqs.(@]), (0), and (|7j), one can readily 
get the desired spinor products in Q, like 

u{p2, A2) ® u{pi, Ai) = u{ko, -X2)u{ko, -Ai 



a/2P2 ■ ko ' ' y/2pi ■ ko 

'2 ± rn2 



:{(l + AiA2)-(Ai + A2)75 



+ ^i[(Ai - X2) - (1 - AiA2)75]} ^0 /;,^""^ . (8) 

2V2pi ■ ko 

In fact (jHl) was obtained in [TTl , but it was not realized that (jS)) can be directly applied 
to the calculation of helicity amplitude, rather, there it was used for evaluating the 
segments of a splitted fermion string. This finding makes the computation of practical 
processes much simplified from the method of Ref. 



3 Applications and Concluding Remarks 

In arriving at eq.® we have taken no special provisos on the fermion mass, hence, it 
will remain valid for massless case. From (jH} we can easily obtain the known identities 
of the inner product of two massless spinors like 

uip, -)u{p', +) = {u{p', +)u{p, -)]* , (9) 



u{p, +)u{p', -) = -u{p', +)u{p, -) , (10) 

and 

\u{p,+)uip',-)\'' = 2p-p' . (11) 

When taking a specific choice for auxiliary vector ko, that is 

^0 = ^ip' - IpI) , ko = 4(|p| - p°)p , (12) 

where a is an arbitrary parameter and p the unit vector of space components of mo- 
mentum p, the polarization vector now becomes the conventional helicity vector: 

n={ — , — p . 13 

m m 

Then, from (jH)) one easily finds that at helicity basis, the spinors u{p, A) satisfies the 
usual projection relation 

u{p, XMp, A) = (1 + fi) ^ ^ ^^4^ 



A 



As the last example, we compute the helicity amplitude for fermion line with a single 
vector current insertion. That is 

A^{p, A; p', A') = u{p, X)i^u{p', A') . (15) 

This kind of amplitude appears quite often in fermion production and decay processes, 
e.g., e"*" e~ //. In the the center-of-mass(CM) system, with choosing 

p = {po, 0, 0, ps) , p' = {po, 0, 0, -ps) , ko = (1, 0, 0, -1) , h = (0, 1, 0, 0) (16) 

and using © the helicity amplitude can be can be evaluated as 

A^{p, A; p', A') = ^1 ^ (s 6^0 - (s - Am^) 6^^ - sko^) 

+ (A' - X)mVs - 4m2 -1(1 - AA')mv^ 6^2} ■ (17) 
Here, m is the fermion mass, s is the CM energy square. From ()17|) we have 

A^{p, \-p\ A') = 2m for A = A' , (18) 



A^{p, A; p', A') = -AVs - 4m2 + Iv^ 6^2 for A = -A' , (19) 

which are exactly the same as given in eq. (23) of Ref. after taking a specific 
momentum direction there. 

In summary, we provide a new approach for evaluating Feynman diagrams symbol- 
ically, which is a straightforward and efficient promotion of the methods developed in 
Refs. |im illj. We have presented examples in displaying the simpleness of this approach, 
and the new method is applicable to both massive and massless fermion amplitude cal- 
culations. Especially, since there is no restriction on the working frame, or in other 
words the fermion momentum directions, by this approach, it is quite convenient to 
encode the helicity amplitudes into a symbolic manipulating program for computer. 
Another advantage of the method is that the expression (jHJ is Lorentz index covariant, 
although one needs to specify the auxiliary vectors in the end for one's convenience, 
which makes the intermediate steps of calculation not so tedious. In the end, we would 
like to mention that most phenomenological applications of massive fermion amplitudes 
can be improved by also including the decay of the fermions. It is often possible to fac- 
torize the problem into production and decay amplitudes, and the new technique could 
be used for both in principle. However, to do this much care should be paid to the 
helicity state conventions. 
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